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1 Introduction 

The origin of CP violation in nature is a very important issue in particle physics since CP 
violation is indispensable for the baryon asymmetry in universe. In the standard model 
(SM), there are two CP violating parameters; the phase of the CKM matrix and the QCD 
theta (6) parameter. The former induces CP violation in flavor changing processes, such as 
CP violating K and B meson decay modes whereas the latter induces flavor conserving CP 
violation, such as neutron electric dipole moment (EDM). The experimental upperbound 
on the neutron EDM gives a strong constraint on 9, as |^| ^ 10^*^^°^"'^^-'. On the other 
hand, the recent measurements of CP asymmetry in the B decay modes at the Babar and 
Belle experiments confirm that the phase of the CKM matrix is the dominant source of 
the CP violation in the K and B decays. It is known that the phase of the CKM matrix 
is not much enough to explain the baryon asymmetry in universe. Therefore, it is very 
important to search for new CP violating phases. 

The minimal supersymmetric SM (MSSM) is one of the most attractive models which 
describe physics beyond the SM. In the MSSM many new CP phases may be introduced 
in both the fiavor-diagonal and changing soft SUSY breaking terms and these new phases 
can contribute to the CP violation at low energy experiments. Recently, the Belle col- 
laboration announced that CP asymmetry in B ^ (pKs is —0.96 ± 0.50, which is 3.5 a 
deviation from the SM prediction [1]. On the other hand, the Babar result on the CP 
asymmetry in B ^ (pKg is consistent with the SM prediction [2]. Many papers appear to 
explain the anomaly in SUSY models [3, 4]. One of the most interesting possibilities is 
that the deviation is a signal of CP violation in right-handed bottom and strange squarks 
mixing. In particular, in the SUSY SU(5) GUT with right-handed neutrinos, the large 
right-handed bottom and strange squarks mixing may be induced due to the large neu- 
trino mixing [5]. However, we have pointed out that the CP violation in right-handed 
squark mixing is strongly constrained by the ^^^Hg atomic EDM through the chromoelec- 
tric dipole moment (CEDM) of the strange quark [6]. 

Motivated by the observation that strange quark contribution to ^^^Hg atomic EDM 
may give an important implication to other observables, we revisit various hadronic EDMs 
induced by the quark CEDMs and the QCD theta term, especially paying an attention 



to the strange quark CEDM. There are several theoretical approaches to calculate the 
hadronic EDMs induced by the quark CEDMs so far. In particular, the hadronic EDMs, 
including the neutron EDM, are evaluated in detail from the QCD sum rule [7, 8, 9, 
10, 11, 12]. In this paper, we take a rather "traditional" approach based on the chiral 
Lagrangian. While the QCD sum rule is a good tool to dictate the QCD dynamics, it is 
difficult to incorporate the sea quark dynamics, especially in the neutron EDM. In fact, 
the sigma term in the chiral perturbation theory suggests that the sea quark dynamics is 
important in the baryon physics [13]. Thus, we adopt the chiral Lagrangian approach in 
this paper, and incorporate the strange quarks in it. We derive the CP violating effective 
Lagrangian with SU(3) chiral Lagrangian in the presence of the QCD theta term and the 
quark CEDMs for the evaluation. 

From this Lagrangian we find that the strange quark CEDM contribution to the ^^^Hg 
EDM is induced by the rj^-ir^ mixing originated from the isospin violation. While the 
strange quark CEDM contribution is evaluated by the eta exchange contribution in the 
previous paper [15], the recent calculation for the Shiff moment [14] suggests that the 
contribution is suppressed by 0(10^^). However, we show that the derived constraint 
on the strange quark CEDM is comparable to the previous result from the eta exchange 
contribution. Also, for the neutron EDM, the loop calculation shows that the constraint 
on the strange quark CEDM is stronger than that from ^^^Hg EDM. The improvement 
of the deuteron EDM is proposed recently [16], and the sensitivity may reach to do ~ 
(1 — 3) X 10~^^ecm. We discuss the sensitivity to the CP violating parameters in QCD, 
including the strange quark CEDM. 

We also discuss the SUSY contributions to the hadronic EDMs through the quark 
CEDMs. We consider the gluino contribution with the mass insertion approximations 
and show the constraints on the CP phases of the flavor changing SUSY breaking terms, 
which are relevant to the SUSY flavor physics, from the hadronic EDM experiments. 

This paper is organized as follows: In Section 2 the CP violation at quark level up 
to the dimension flve operators is reviewed. In Section 3 we derive the effective nucleon 
interactions induced by the CP violating interaction at quark level. In Section 4 we esti- 
mate various hadronic EDMs with the CP violating nucleon interactions. In Section 5 we 
consider SUSY contributions to the hadronic EDMs with the mass insertion approxima- 



tion and show the constraints on the SUSY parameters. Also we discuss the correlation 
between the hadronic EDM and the CP asymmetry in B ^ <pKs- Finally, Section 6 is 
devoted to conclusion and discussion. 

2 Hadronic CP violation in quark level 

The CP violation in the strong interaction of the light quarks is dictated by the following 
effective operators, 

l^ = ^^GG^ Y: tfqgsiGah.q, (1) 

q=u,d,s 

up to the dimension five ones. Here, G^,^ is the SU(3) gauge field strength, G^,^ = 
lef.^p^f"' and Ga = G^^^a'^'T''. The first term in Eq. (1) is the effective QCD theta 
term, and the second term is for the quark CEDMs. The 9 parameter can be 0(1) 
generically, however, it is strongly constrained by the neutron EDM experiments. One 
of the most elegant solution is to introduce the Peccei-Quinn (PQ) symmetry [17], since 
the axion makes 6 vanish dynamically. On the other hand, when the quark CEDMs are 
non-vanishing, the theta term is induced again even if the PQ symmetry is introduced, 
as pointed out by Ref. [18]. Thus, in the evaluation of the effect of the quark CEDMs, 
we need to include the QCD theta term systematically. In this section we review the 
CP violating effective interactions in the strong interaction of the light quarks up to the 
dimension five operators and the role of the PQ symmetry. 

We start from the case without the PQ symmetry, first. Let us take a basis, in which 
the QCD theta term is vanishing, by the chiral rotation as 

£<^ = - J2 rngqiaq-f5q+ J2 ^^q9s{Ga)-i5q, (2) 

q=u,d,s q=u,d,s 

where we have used |^| <^ 1. The first term is induced by the axial anomaly and a^'s 
satisfy the following relations, 

^= E «r (3) 

q=u,d,s 



Furthermore, we also impose for convenience a condition that the CP violating tadpoles 
for 7r° and if should vanish, (tt^ \Cgp\ 0) = {if \Cg^\ 0) = 0. Using the PCAC relation, 
these conditions and Eq. (3) determine a^'s as 
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where itIq is the ratio of the quark-gluon condensate to the quark one. From the QCD 
sum rule, it is estimated as [19] 
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Notice that the CEDM of the strange quark gives sizable contributions to au and ad when 
du/rriu ~ dd/rrid ~ ds/iris. 

This situation is changed when the PQ symmetry is introduced, since 9 is promoted 
to a dynamical field, axion (a), 

C = a^GG. (8) 

87r 

If the quark CEDMs are vanishing, ^(= (a)) is aligned dynamically to the zero at the 
minimum of the axion potential. However, the axion potential is modified in the presence 
of the quark CEDMs [18], 



Kflfla) = Kia+ i^Ka^, 



where 



K 






limi / rf^a;e^^^(0 



rv ~ (y ~ 

—GG(x) —GG(0) 



0), 



ry ^ , r] 

-^GG{x) Y: ^^ Q 9s{Gah, q{0) 
on , I 

q=u,d,s 



(9) 

(10) 
0). (11) 



The linear term with respective to a, which is proportional to the quark CEDMs, induces 
non- vanishing vacuum expectation value for a. The matrix elements, K and Ki, are 



evaluated using the current algebra [20, 18], and they are given as 

K = . """'^ . (0|gg|0), (12) 

Ki = -7^7^^^^^^ E ^{0\q^9siGah,q\0). (13) 
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q=u,d,s y 

Thus, the minimum of the axion potential is shifted in the presence of the quark CEDMs 
and the following 6 parameter is effectively induced as 

0.J=/a)) = -^ = ml y ^. (14) 

q=u,d,s y 

Plugging Eq. (14) into Eqs. (4-6), it is found that a^'s become simpler under the assump- 
tion of the PQ symmetry, 

tty = ^0^ {q = u,d,s). (15) 

The PQ symmetry suppresses the contribution of the strange quark CEDM to «„ and a^,- 
This means that the valence quarks themselves in nucleon do not suffer from the strange 
quark CEDM. 

3 CP violating nucleon interactions 

In the previous section we have considered the CP violating interactions at the quark level. 
We need to translate those quark level interactions into hadronic interactions in order to 
calculate the hadronic EDMs. This is a rather difficult task, and the sizable theoretical 
uncertainties are expected for evaluation of hadronic matrix elements. Here, we consider 
an approach based on the chiral Lagrangian technique. It is possible to discuss hadronic 
CP violation with the effective Lagrangian in a systematical and simple way. In this 
section we derive the effective CP violating nucleon interactions using the SU(3) chiral 
Lagrangian in order to incorporate the strange quark contribution in nuclear interactions. 
First, let us review the chiral Lagrangian without CP violation. We consider the 
SU(3)LxSU(3)i? chiral symmetry and the effective Lagrangian is written in terms of the 



meson fields 
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The effective Lagrangian invariant under SU(3)LxSU(3)i{ is given by 
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with U = ^'^. The phenomenological parameters D and F are 0.81 and 0.44, respectively. 
The effective Lagrangian is expanded by the derivatives, and the lowest order terms are 
important at the low energy. The terms relevant in our paper are summarized in Appendix 
A. 

The pseudoscalars are massless due to the SU(3)LxSU(3)ii: symmetry. However, the 
SU(3)lxSU(3)r symmetry is violated by small quark masses and the pseudoscalars be- 
come massive. Also, the isospin breaking effect has an important role in the EDMs of the 
^^^Hg atom and deuteron as shown later, since it leads to the n^-i]^ mixing. 
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where tt' and rj' are for the mass eigenstates. The mixing angle is given as tan^ 
-v^(m„ - md)/(4ms) ~ 0.01 for m„ = 5.1 MeV, ma = 9.3 MeV and m^ = 175 MeV. 



Now we introduce the CP violating terms to the chiral Lagrangian. We need to eval- 
uate the matrix element {BaTv'^ KAtpI ^b) in order to obtain the CP violating interactions 
between nucleons and pseudoscalars induced by the CP violating terms in Eq. (2). It is 
reduced by the PCAC relation, 

{BaiT''\i/:gp\Bb) = Yl -^[{,aqmq + aq>mq,){Ba\qTcq\Bi,) 



q,q'=u,d,s J ''^ 

(21) 



:^ + ^){BaW9s{Ga)TMB,) 



where Tc is a generator for the flavor SU(3) symmetry. 

The matrix elements of the scalar operators in the right-handed side in Eq. (21) are 
represented as 



{Bamq\B,) = XTi [BaB,T,)+YTT[BaT,B,), (22) 

{Ba\qlq\B,) = S^i^Z, (23) 

from consideration of the transformation property under SU(3)LxSU(3)jt:. Here, 1 is the 

(3 X 3) unit matrix. The phenomenological parameters X, Y and Z can be determined 

by the baryon mass splittings and the nucleon sigma term [13]. 

1 

X = (-2m= + 3mA - m^) = -idd) + (^s), (24) 

2ms 

Y = — (m= - ms) = -(dd) + (uu), (25) 

rus 

^ 3(7 3 m= — vriK 

^ = ;;^:t^- 2^^ =<''''>+<''''> +<'^>' '2«' 

where (gg) = (p|gg|p) and o represents the nucleon sigma term, which is estimated as 
o" ~ 45 MeV. Numerically they are given by 

X = -1.35, Y = 0.72, Z = 7.7, (27) 

{uu) = 3.5, (dd) = 2.8, {ss) = 1.4. (28) 

For the evaluation of {B^ \q'gs{Go')Tcq\ Bf,), we use the following relation adopted in 
Ref. [15], which is inspired by the QCD sum rule, 

{Ba \q9siGa)q\ B,) ~ ^m^(5, |gg| B,) (29) 



for each quark. From Eqs. (22) and (29) 

{Ba\q9s{Ga)T,q\Bh) = ^ml [XTi (KB.T,) +YTt [KT.B,)) . (30) 

This relation is, off course, approximate one, however, it makes the formula simple as 
shown now. 

From the above consideration, Eq. (21) can be reduced as 

{Ban'\tjC<jp\B,) = ■^{Ba\q{T,,A}q\B,). (31) 

Jit 

The matrix A in Eq. (31) is a diagonal matrix, A = diag{Au, A^, Ag), with the compo- 
nents 

Ag = agnig-^-ml (32) 

Therefore, the CP violating Lagrangian is obtained in terms of the baryon and meson 
fields, 

Atp = -^{xTt[bB{M,A})+YTt[b{M,A}B 

+ -{Z -X -Y)Ti{AM)Ti(bb)Y (33) 

The relevant terms to our discussion in Eq. (33) are summarized in Appendix B. 

The matrix elements Ag^s are written in terms of the quark CEDMs and the QCD 
theta parameter. When there is no PQ symmetry, they are given by 

A^ = 3.3 X lO"-'^^ - 0.53du - O.Udd - 7.5 x 10"^4 GeV, (34) 

Ad = 3.3 X lO"-'^^ - 0.26du - 0.41 J^ - 7.5 x 10~^4 GeV, (35) 

As = 3.3 X lO"-'^^ - 0.26rf„ - O.Udd - 0.274 GeV. (36) 

Here, rfg's are written in GeV unit. On the other hand, when the PQ symmetry is 
introduced, Ag^s are written in terms of the quark CEDMs and become rather simple, 

Ag = -0.27dg GeV, {q = u,d,s). (37) 



4 Hadronic EDMs 

In this section we consider the hadronic CP violation. Among various observables, the 
EDMs of neutron or atoms are very sensitive to the flavor conserving CP violation. In 
particular, the experimental upper bound on the EDMs of neutron [21] and ^^^Hg atom 

[22], 

I^Hgl < 2.1 X lO^^^e cm, (38) 

|rf„| < 6.3 X 10~26e cm, (39) 

give strong constraints on hadronic CP violation. Also, the improvement of the deuteron 
EDM is proposed recently, and the sensitivity may reach to rf/j ~ (1 — 3) x lO'^'^e cm [16]. 
Therefore, we estimate the EDMs of ^^^Hg atom, neutron and deuteron with the chiral 
Lagrangian obtained in the previous section. 

4.1 i^sjjg Atomic EDM 

The ^^^Hg atom is a diamagnetic atom, in which electrons make a close shell. In such 
atoms, the atomic EDMs are primary sensitive to the CP violation in nucleons and rep- 
resented by the Shiff moments (S). The Shiff moment generates the T odd electrostatic 
potential \4ff = AirS (/ ■ V)6{r) with the nucleus of spin (/). The ^^^Hg atomic EDM is 
calculated in terms of the Shiff moment S [23] 

c^Hg = -2.8 X 10"^^ 5(^^^Hg) e cm, (40) 

where S is in unit of e /m^. This implies that the Shiff moment is bounded as [24] 

\S\ < 0.75 X 10-" e fm\ (41) 

from Eq. (38). 

The Shiff momentum is induced by the pion and eta exchanges with the CP violat- 
ing coupling. From the formula derived in the previous section, the CP violating pion 
interactions are given as 

^^^ = ^({A^ + A,)-^{A^-A,)\{{uu)-(dd))NT^N7,^ 



+ -4^ (A„ - A,) ({uu) - (dd)) (Nr^Nn- - SNt^Ntt^) 
6v3/„ ^ ^ 

= g^^l^NT'^NTT'^ + g^^LNNir' + g^^l^ {Nt^Nt,- - SiVr^iVTr") , (42) 

where A^ = {p,nf , and 'qIIin, 9-kNN ^"^^ 'SInn correspond to the isoscalar, isovector and 
isotensor couphng constants, respectively. Here, we include the tt^-i]^ mixing. The CP 
violating eta ones are 

£5^" = -^(^{A^ + A,){{uu) + (dd))-2A,{ss)^NNv (43) 

^2717 ^^" ~ ^""^ (^^""^ ~ ^^"^^^ '^^^^V 
^ g^^l^NNr] + g%NnNr]. (44) 

The recent evaluation for the Shiff moment of ^^^Hg [14], in which the many-body 
treatment is performed, reveals that the core polarization effect reduces the isoscalar and 
tensor channel contributions while the isovector channel one is comparable to the previous 
results [25]. Their result for the Shiff momentum is^ 

-^ ^ X 10-2 / (0) (1) _ ^ (0) (1) X 

O.O X iU X \g.„NN9TiNN 2^vNN9riNN) 

+9 X 10-3 X ^S^^S^v [e fm% (45) 

Here, glpfj^ and gliNj^ is the CP even pion and eta coupling constants, respectively,^ 

CcP = '^{D + F)-N^,T'^Nn'^'-^{D-3F)-N^,Nr]'' (46) 

^ g^^l^N^,T^N7r^'+g%N^,Nr^''. (47) 

The results in Ref. [25], which neglect the core polarization, correspond to S'(^^^Hg) = 

-0.086 X g'^ifNig^'^NN + gnNN - '^9wNn) e fm^ for the pion contribution [14]. 

^Whilc only pion exchange is included in in their calculation, they show that the correction from 
the finite range of the pion interaction is small. Thus, the above introduction of the eta contribution is 
justified. 

^The isovector coupling in the CP even pion interaction is also induced by the ?7°-7r° mixing, but we 
neglect it, since it does not change our qualitative discussion. 
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From this evaluation, it is found that the EDM of ^^^Hg atom is sensitive to (A„ — Ad) 
through the isovector channel, and the contribution proportional to (Au+Aa) is suppressed 
by 0(10^^) due to the rj^-n^ mixing or the isoscalar channel. This means that the ^^^Hg 
atomic EDM is insensitive to the QCD theta parameter. Also, notice that 'gUfj^^ depends 
on As through 71^-1]^ mixing. While eta itself has the CP violating interaction which 
depends on Ag, the eta exchange contributions are suppressed due to the isoscalar channel 
or the mass. Thus we concentrate on the isovector channel of the pion exchange from the 
CP violating pion interaction with non-vanishing rj^-TC^ mixing. 

The Shiff moment from non- vanishing yl^lfj^ is represented in terms of Ag's as 

^(i99Hg) = -(23 X (A„ - Ad) + 0.13 X A,,) e fm\ (48) 

where Aq{q = u, d, s) are given in GeV unit. The current experimental bound on the 
EDM of ^^^Hg atom implies that 

e\du -dd + 0.005141 < 2.4 x lO'^^e cm. (49) 

This is almost independent of whether we impose the PQ symmetry or not, since it 
depends on A^ — A^. The strange quark CEDM is bounded as 

e\ds\ < 4.7 X 10~^^e cm. (50) 

It is found that this bound does not have a qualitative difference, especially in the 
contribution of the strange quark CEDM, from one based on Ref. [25], that is, e\du — dd + 
0.012(is| < 7 X lO^^^e cm [15]. In our evaluation the contribution of the strange quark 
CEDM comes from the vj^-ti^ mixing, and it is suppressed by ~ (m„ — md)/{msm\). On 
the other hand, the eta exchange contribution in the previous evaluation is suppressed by 
~ l/""^^- Thus, it is found that they are comparable. 

4.2 Neutron EDM 

The evaluation of the neutron EDM has a long history. However, the precious calculation 
is still a very difficult issue. One of the simple estimations is based on the Naive Dimen- 
sional Analysis [26] . However, we can know only the order of magnitude of the neutron 
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EDM in this method. More elaborated calculations are done based on the QCD sum rule 
[7, 8, 9, 10, 11, 12]. The QCD sum rule is an established technique, and it has a merit for 
the evaluation of the neutron EDM. In Ref. [11] the QCD sum rule analysis reproduces 
the ratio of the quark EDM contributions in the neutron EDM, which is expected from 
the non-relativistic quark model. On the other hand, they argue that the PQ symmetry 
suppresses the strange quark CEDM contribution in the neutron EDM. However, the sup- 
pression is for the valence quarks, not for the sea quarks. They introduce the contribution 
of the strange quark CEDM only through the vacuum expectation value of qi'^^q with q 
the valence quark under the CP violating background. It should be suppressed under the 
PQ symmetry as discussed in Section 2. In fact, the sigma term in the chiral perturbation 
theory suggests that the sea quark dynamics is important in the baryon physics [13]. As 
shown in Eq. (28), the matrix element of the strange quark is comparable to those of the 
other light quarks in nucleon. 

Thus, we evaluate the neutron EDM from the effective Lagrangian obtained in the 
previous section by the traditional loop calculation so that we include the strange quark 
contribution. The neutron EDM ((i„) is defined by 

-Cedm = — ^ncr^^-f^nF^". (51) 

The EDM is induced by the one loop diagrams with the charged mesons in Fig. 1 and 
they are written as 

dn = j^, {A^x^:^ + A,x^-^ + Ax^:^^) , (52) 



47r2/.^ 



where 



(n) _ (D-L P\ (/^oA - /rlrl\\]^^!!]2. J- (D - P') f/rlrl\ ^ /^o\\u^I!]Z. 



^U 






{D + F)( {uu) - (dd) )\og^ + {D-F)( {dd) - (ss) log ^^ = 2.0, (53) 
^ ^ m^ ^ ^ rriK 

(D + F) ((uu) - (dd)) log ^ = 1.7, (54) 



a;(«) = {D-F)[{dd)-{ss)) log ^^ = 0.33. (55) 

Here, we have considered the most singular terms, which have logarithmic singularities in 
the chiral limit originated from n and K mesons. 
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Let us consider the constraint from the neutron EDM experiment. When there is no 
PQ symmetry, the neutron EDM is given in terms of the quark CEDMs and the QCD 
theta parameter, 

rf„ = (7.7 X 10"^¥ - 4.6 X J„ - 3.0 x J^ + 0.33 x 4) e cm. (56) 

If we impose the PQ symmetry, 

dn = (1.6 X du+ 1.3 X dd + 0.26 x ds) e cm. (57) 

Here, the quark CEDMs in the above equations are written in cm, unit. For comparison 
to the result of the QCD sum rule in Ref. [11], the contributions from the up and down 
CEDMs are slightly larger in our evaluation. It is reasonable since the sea quarks also 
contribute to the neutron EDM in addition to the valence quarks. 

From the current experimental bound of the neutron EDM in Eq. (39), we obtain the 
following bounds for the the theta parameter and the quark CEDMs for a case that the 
PQ symmetry does not exist, 

1^1 < 8.2 X 10"^\ 

e\du\ < 1.4 X 10"^^ e cm, 

e\dd\ < 2.1 X 10~2^ e cm, 

e|4| < 1.9 X 10"^^ e cm. (58) 

If the PQ symmetry works, 

e\du\ < 3.9 X 10^^^ e cm, 
e\dd\ < 4.8 X 10"^^ e cm, 
e\ds\ < 2.4 X 10"^^ e cm. (59) 

Here, we assume no accidental cancellation among the various contributions. 

Notice that the constraint on the strange quark CEDM is one order of magnitude 
stronger than that from the ^^^Hg atomic EDM, though the prediction for the neutron 
EDM may be expected to have more theoretical uncertainty. While the experimental 
bound on the neutron EDM is weaker than that of the ^^^Hg atomic EDM, the strange 
quark constraint is not suppressed by the rj^ mass or the isospin violation. 
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4.3 Deuteron EDM 

The new measurement for the deuteron EDM is proposed in Ref. [16], and it is argued 
that one or two orders of magnitude improvement may be achieved relative to the current 
bounds on the QCD theta parameter and the nuclear force. It may give a stringent 
constraint on the new physics or discover the signal. Also, the nuclear dynamics in 
deuteron is rather transparent and the theoretical uncertainty is expected to be small. 
The deuteron EDM is given by 

dD = {dn + dp) + d^"" , (60) 

where dp is the proton EDM and the second term comes from the CP violating nuclear 
force. 

In Ref. [27] they show that the chiral logarithms in the sum of the proton and neutron 
EDM exactly cancel in the chiral perturbation theory, and they adopt the QCD sum rule 
for it. However, the cancellation in the chiral perturbation theory does not come from 
some symmetry. The chiral logarithms in the neutron and proton EDMs come from loop 
diagrams of charged mesons in which photons are attached to the changed mesons. Thus, 
it is obvious that the chiral logarithms are canceled in the SU(2) chiral perturbation 
theory. However, it is not true if we introduce the strange quark and the corresponding 
mesons. 

In the SU(3) chiral perturbation theory, the sum of the proton and neutron EDMs is 

(dp + dn) = -47^(A„4"+P)+A,xi"+^)+^x("+^)), (61) 



where 



47r2/^ 



>+^) = --( D({uu)-5{dd)+4{ss)) + 3F((uu) + (dd)-2(ss))) log ^ 
3 ^ ^ rriK 

+-D{6{uu) - 5{dd) + (ss)) + F{4{uu) - (dd) - l^s)), (62) 

i"+^^ = iz^(3(llM)-10(rfrf) + 3(ss)) + F((llM)-(ss)) (63) 

= — (DiiJlu) - hCdd) + 4(ss)) + 3F((mm) + Cdd) - 2(ss))) log ^ 
3 ^ ^ rriK 

+-D{{uu) - 5(dd) + 6{ss)) + F{{uu) + (dd) - 4(ss)) (64) 



^d 
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and then x^"'^*'^ = 6.4, x^ = —2.7, and a;^"+*') = —0.40. The logarithmic terms in the 
proton EDM come from the K^-A'^ loop in addition to the K~^-J]~ and 7r+-n loops. We 
also include constant terms which are not suppressed by proton mass. The strangeness 
has an important role in the deuteron EDM. When the PQ symmetry does not exist, 

{dp + dn) = (6.4 X 10"^^^ - 7.7 x J„ + 0.73 xdd + 0.23 x J,) e cm. (65) 

If it does, 

(dp + dn) = {— 5.1 X du + '2.1 X dd + 0.32 X ds) e cm. (66) 



In Ref. [28] it is shown that d^^ comes from the isovector coupling of pion gl^^pf, 
given in Eq. (42), as 

<^ = -2.1 X W' X g^g'^^l^ e cm. (67) 

It implies 

^nNN ^ (-11 X du + 11 X dd- 0.063 xds)e cm. (68) 

From the above evaluation, the measurement of the deuteron EDM will be a stringent 
test for the SM. If they establish the sensitivity of do ~ (1 — 3) x lO^^'^ecm, they 
can probe ed^ ~ edd ~ 10"^^ e cm, from the nuclear force as shown in Ref. [27], and 
9 ~ 10~^^ e cm, from the nucleon EDMs. Also, for the strange quark CEDM, they may 
reach to eds ~ 10^^^ e cm from the nucleon EDMs. 

5 SUSY contributions 
5.1 Hadronic EDMs 

Let us consider the SUSY contributions to the hadronic EDMs. In the MSSM many 
CP violating parameters can be introduced in the soft SUSY breaking terms. If there are 
SUSY CP phases, the quark CEDMs are induced through one loop diagrams. The bounds 
on the SUSY contributions to the hadronic EDMs can be estimated from the results in 
Section 4. 
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In our previous paper [6], we have considered the gluino contribution to the strange 
quark CEDM induced by the flavor mixing in the MSSM. While the EDMs are flavor 
conserving phenomena, the stringent bounds for the EDMs give constraints on even the 
flavor mixing, and they have an impact on the flavor physics, such as B meson decay. Thus, 
we consider the gluino contributions including the flavor violation inside the Feynman 
graph in order to demonstrate the constraints on the SUSY CP phases from the our 
calculations.^ 

In general SUSY models, both left-handed and right-handed squarks have flavor mix- 
ings. In this case, the CEDM of the i-th light quark qi is generated by a diagram in 
Fig. 2(a), and it can be enhanced by rnqjiriq. when j > i, {i,j = 1, 2, 3). Using the mass 
insertion technique, it is given by 

4 = c^3 (-^^i(^) - 3iV2(x)) Im [J^f] , (69) 

where rrig and rriq are the gluino and averaged squark masses and c is the QCD correction, 
c ~ 0.9. The functions Ni are given as 

,, , , 3 + 44a;-36a;2-12a;3 + a;^ + 12x(2 + 3a;)loga; 

"^^^"^ = 6(^3Tji ' (^0) 



10 + 9a; - ISx"^ - x^ + 3(1 + 6a; + Sx"^) logx 
3(a;-i: 



^2(a:) = :^Yz T^e • (^1) 



The flavor violation in the squark mass terms contributes to the quark EDMs via combi- 



nations of JjJ , 



Aq) _ fdq)^ fdq)-s ( dq) 



W = ('^i'/).. ('^iU. (4U.- (72) 



(l)\. {X(l)\. ^^A fxM 



The mass insertion parameters, {S^Dij, {^Rjijij and {oj^lijij, are defined as 



dq)\ _ \JiJij ..(g) ^ _ V iR/ij 



i^LLJij — ^2 ' i^RRJiJ ~ ^2 ' \''^) 



m| m| 



for q = u,d, and 



_ m,^.(Af -/xtan/3) m,^ {a^^ - f^ cot p) 

K'^LRJij — "ij 2 ' \"LR)ij — "ij ^2 ' V ' ^/ 



■^In Ref. [29], it is pointed out that the chargino diagrams can also give large contribution to ^^^Hg 
EDM. A comprehensive analysis of the SUSY contributions to the hadronic EDMs will be discussed 
elsewhere. 
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Here, {fn^^ ) is the left-handed (right-handed) squark mass matrix. Here, we assume 
that the flavor-conserving SUSY breaking terms do not have the CP phases since they 
are aheady strongly constrained, and we keep terms induced by the off-diagonal terms in 
the squark mass matrices 

In the typical SUSY models, the left-handed squark mixings are governed by the CKM 
matrix, and even if the universal scalar mass hypothesis is imposed at high energy scale, 
the radiative correction induces the off-diagonal terms as 



(4^)13 = Oi\') ^ 8 X 10-^ 



(4^23 = 0{\') ^ 4 X 10 



1-2 



(75) 



for A ~ 0.2. On the other hand, the right-handed squark mixing is rather model depen- 
dent. In the SU(5) SUSY GUT with right-handed neutrinos, large right-handed squark 
mixings can be induced by the neutrino Yukawa couplings [5] . They are constrained from 
the K^-'K^, D°-W, B^-W mixings experiments [30]; 





< 


4 X 10"^ 


i^RR)l2 


< 


1 X 10"\ 


{^RR)13 


< 


1 X 10^\ 



(76) 



for rriq = rrig = 500 GeV. 

The quark CEDMs are estimated from Eq. (69) as 



€Ctu 



edri 



-7.9 X 10^^* sin ^(^^^e cm x 

(41)l2 



-28, 
-3 



V500GeV/ \ 3 x 10-^ 
5.7 X lO-^Sin^P^ecm x 



ma 



V500GeV/ 
-4.0 X 10^^^ 



V500GeV/ 
- 9.6 X 10~2^ 




0.1 



(4t)3 1 

0.1 



(^Rij)21 



4 X 10-2 




ficot (3^ 



500GeV 



At — ficotp^ 
500GeV 



/ fitanP \ 
1^5000Gevj 



(77) 



sin 9^'' e cm x 
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TTlg 



500GeV/ 



(4t)l3 

8x 10- 



{^rrJsi 



I yutan/3 \ 



0.1 I vsoooGev; 



ed,. 



4.8 X lO'^^sin 

3 



SOOGeV/ 



^(3) 



ecm 



LL)23 

0.04 



[ "Rr)3 2 
0.1 



/ /itaii/S \ 
UoOOGevJ 



(7J 



(79) 



where we take m^ 



771; 



and 9^^' is the phase of the SUSY parameters, 



t.e. 



d^'> 



arg[J,(^^ 



u 



Here, we take the experimental bounds for mass insertion parameters of the right-handed 
squarks in Eq. (76) while the left-handed ones are given by Eq. (75). We have neglected 
the A-terms in d^ and dg since they are subdominant. 

When one of the quark CEDMs saturates the hadronic EDM bounds, we can obtain 
the following bounds on the SUSY CP phases from the ^^^Hg atomic (neutron) EDM 



experiments, 

sin^p) 



sm 



^P) 



< 30(47) 



< 4.2(6.5) X 10 




^0) 



(2) 



sin^ 



sin^f 



60(113) 




['JRr)21 

4 X 10-2 
1 



/ fitanf] \ 



V5000GeVy 

-1 



i^RRhi] fjitan£_y 



0.1 



sm 



^P) 



< 0.98(0.048) 



0.04 



0.1 



V5000GeV; 

/ yUtan/3 \ 
VSOOOGeVJ 



(82) 
(83) 
(84) 



where we take rriq 



rrig = 500GeV. The above bounds on 9^^ and 9^ are looser than 



than the stringent K^-K constraint. On the other hand, the CP phases related to the 
1-3 and 2-3 mixing angles are constrained by the hadronic EDMs significantly. These 
bounds are expected to be improved furthermore by the deuteron EDM measurements. 



^Here we do not consider the quark EDM contributions to the neutron EDM. We find that they are 
subdominant compared with those of the quark CEDMs in our calculation. 
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5.2 Hadronic EDMs and CP asymmetry in 5 ^ (I^Ks 

Let us consider a correlation between ds and S^Ks in the SUSY models. In Ref. [6], we 
have shown that there is a strong correlation between them when both left-handed and 
right-right handed squarks have flavor mixings. In such a case, the dominant contribution 
to S^Ks is supplied by a diagram with the double mass insertion of (5^^)32 and (5^^)33 
(Fig. 2(b)). The contribution of Fig. 2 (b) is given as 



where 



H - -C^^m,MGcr)bL, (85) 



^^ = ^^(^S)33(4'i)32(-^Mi(x)-3M2(a:)), (86) 



up to the QCD correction. Here, 



.., . 1 + 9a; — 9a;^ — x^ + (6x + Qx'^) log a; ,„_, 

Miix) = ^^^-^p , (87) 

,./x 3-3x'^ + il + 4:X + x^)\ogx 

M2{x) = 



a;-l 



In a limit of a; — ^ 1, C^ is reduced to 



^" - 6sS^(^S)33(4'i).. (89) 



Using Eqs. (69) and (89), we find a strong correlation between dg and Cq as 

M1)23C^] (90) 



ds = ^ — Im 

47r2 21 



up to the QCD correction. The coefficient 11/21 in Eq. (90) changes from 1 to 1/3 for 
< a; < 00. 

In Fig. 3, we show the correlation between dg and S(pKs assuming a relation dg = 
-mb/(47r2)(ll/21)Im[(4L)23<^8^] up to the QCD correction. Here, we take (4l)23 = 
—0.04, arg[C|^] = 7r/2 and \C^\ corresponding to 10^^ < | (5^^)32! < 0.5. The matrix 
element of chromomagnetic moment in i? ^ (j)Ks is 



{<lyKs\f^m,i7s,a'^'^T,^Pnb,)G%\B,) = «:^(e^p^)/^mjF+(m|), (91) 
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and K = —1.1 in the heavy-quark effective theory [4]. Since k, may suffer from the large 
hadron uncertainty, we show the resuhs for k, = —1 and —2. From this figure, the 
deviation of S^i^k^ from the SM prediction due to the gluon penguin contribution should 
be suppressed when the constraints on dg from the ^^^Hg atomic and the neutron EDMs 
are applied. Comparing Fig. 2 in Ref. [6], the ^^^Hg atomic EDM bound allows a sizable 
deviation of 5*^X8, especially in cases oi k = —2 and |(5)^/)23| smaller than 0.04. It comes 
from the new theoretical estimation of the strange CEDM constraint in Eq. (50) and the 
numerical error in previous calculation. However, we find that the neutron EDM gives a 
strong bound on SfpXs- Moreover, S^Ks may be constrained further by the future deuteron 
EDM measurements. Therefore, the hadronic EDMs give a very important implication 

to S^Ks- 

6 Conclusion and discussion 

We have considered hadronic CP violation induced by chromoelectric dipole moments of 
light quarks and the QCD theta parameter, especially paying an attention to the strange 
quark CEDM. First, we have derived the effective CP violating nucleon interactions in- 
duced by the CEDMs and the theta parameter using the chiral Lagrangian technique. In 
order to take into account the strange quark contributions, we have used the SU(3) chiral 
Lagrangian. Using the effective CP violating Lagrangian, we calculated the EDMs of the 
^^^Hg atom, neutron and deuteron. 

The ^^^Hg atomic EDM is sensitive to the CP violating nuclear force induced by the vr 
and 77'' exchange diagrams. Though the contribution of strange quark EDM is evaluated 
from the eta exchange diagrams in the previous papers, it is found in the new evaluation 
of the Shiff momentum that the isoscalar channel contribution in the ir and rj^ exchange is 
suppressed. We found that the isospin breaking nucleon, interactions originated from the 
^o_^o jYiixing, leads to the similar constraint on the strange quark CEDM to the previous 
one. 

We evaluated the strange quark CEDM contribution to the neutron EDM using the 
standard meson loop calculation. This is originated from the one loop diagram involving 
the K meson. We have found that the neutron EDM gives a stronger bound on the strange 
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quark CEDM in this calculation than the current ^^^Hg atomic EDM experiment, since 
the contribution to the neutron EDM is suppressed by the loop factor at most. While this 
calculation has theoretical uncertainties, it also suggests that the strange quark CEDM 
should be small. 

The new technique for the measurement of the deuteron EDM has a great impact 
on the strange quark CEDM if it is realized. If they establish the sensitivity of d^ ~ 
10~^^e cm, we may probe the new physics to the level of eds ~ 10~^^ e cm, which is 
stronger than the bound from the neutron EDM. 

In order to demonstrate an implication of our result on the SUSY models, we calculate 
the gluino contributions to the quark CEDMs with the flavor violating mass insertions. 
It is usually considered that the EDMs are sensitive to the flavor diagonal CP phases. 
However, when both left-handed and right-handed quark mixing exist, the CEDMs can 
be enhanced by the left-right squark mixings. Since the typical SUSY models have the 
left-handed squark mixing, the EDMs can give strong constraints on the flavor dependent 
SUSY phases. These constraints on the SUSY phases can give important implications to 
other SUSY phenomenology, including the B physics. As an example, we have show that 
there is a strong correlation between the strange quark CEDM and 8,^X3 ■ The current 
bounds on the strange quark CEDM from ^^^Hg atomic and the neutron EDMs imply 
that the deviation of S^^Ks from the SM should be strongly suppressed. 
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A CP even terms in chiral Lagrangian 

In the chiral Lagrangian the interaction can be expanded by derivatives and the low- 
est order terms are important at the low energy. Here, we summarize the relevant CP 
even terms in the chiral Lagrangian. The relevant CP conserving couplings between the 
nucleons and the pseudoscalars are 

= ^{D + F)in,p7r' - ^{D + F)n75n7r° + !!h^±^(D + F)n^,p7r- 
^J^{D - 3F)^5pri^ - -^^{D - 3F)n-f5nri^ 



~ o 7?^ (-P + 3F)A"75nir" + — — {D - F)T.^-f5pK 

- -^-7^^{D + 3F)A0-f5pK +— /|7 — {D-F)J: -i^nK , (92) 

where we have used the equation of motion for the nucleon fields, N'-j^'-j^^N'O^M = 
—i{inj^ + mi\f!)N'~ft^N'M. The coupling constants are given by numerically, 

S'ttpp ^ 12.0, Qnnn ^ 12.0, (j'vrpn ^ il-i), W^) 

gripp = 2.98, grinn = 2.98, gKj^n = 5.98, (94) 

for D = 0.81 and F = 0.44. 

The pseudoscalars are massless due to the SU(3)LxSU(3)ij symmetry. However, the 
SU(3)i,xSU(3)ij symmetry is violated by small quark masses and the pseudoscalars be- 
come massive. When the quark masses m = diag{mu,m(i,ms) are taken into account, we 
can introduce the following terms, 

£i = v^Ti (u^m + mU 

+ aiTr (b [^^m^^ + ^m^) b) + asTr (bB {C^m^^ + ^m^ 
+ feiTr (^75 U^m^^ - ^m^) b) + 62 Tr (575^ U^m^^ - ^m^ 



-2^Tr (^M^rnj + 2aiTr (^mfi) + 2a2Tr (^BBm 

J n 

--L^Tr (B-f5{Mm + mM)B) - ^^Tr (B-f5B{Mm + mM)) + ... . (95) 
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In the estimation of the hadronic EDMs, the isospin breaking effect is important as shown 
in text. The isospin symmetry is violated by the quark mass term and the tt^-i]^ mixing 
occur through the mass terms. From Eq. (95), the 71^-1]^ mass matrix is given by 

^-^'" /2 ^""'"^^ [ -^^im^-md) |(m„ + m, + 4m,) )[v )' ^ ^ 

The mass eigenstates are defined by 

\ri' J \ -sin^ cos^ J \ V J ' 

where tan^ = —^/3{mu — ma) / {4:ms) ~ 0.01 for m„ = 5.1 MeV, nid = 9.3 MeV, rus = 175 
MeV. 

B CP-odd terms in chiral Lagrangian 

From the effective Lagrangian in Eq. (33), the CP violating nucleon interaction terms are 
written as follows. Here, we do not include the tt^-i]^ mixing. 

= — (^Au{uu) - Ad(dd)j ppiT^ 
+ — (Au{dd) — Ad{uu)) nriTT^ 



u 

+ —f=-r {A^ + Ad) Uuu) - (dd)) npir' 
+ ^=r- (Au {uu) + Ad (dd) - 2 As {ss)) ppr]'^ 
+ -=r- ( A„ (dd) + Ad {uu) -2As{-ss)) nnrf^ 
^ {Ad + As)((dd)-{ss))W^pK^ 



V2U 
+ -L^Ad + As)((dd)-{ss))TPnK^ 

+ -^ {Ad + As) ({dd) + {ss) - 2{uu)) l^nW. 

2VojTr 

- -L(^A^ + As){{dd)-{ss))TPpK- 
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+ 



1 



2V3U 
1 



{Au + As) Udd) + (ss) - 2{uu)) A<^pK 



V2U 



(A„ + As) (dd) - (ss) T.-nK~. 



(9^ 
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Figure 1: One loop diagrams for the neutron EDM. The blob represents the CP violating 
coupling. 
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Figure 2: (a) The dominant diagram contributing to the CEDMs of light quarks when 
both the left-handed and right-handed squarks have flavor mixings, (b) The dominant 
SUSY diagram contributing to the CP asymmetry in B ^ (pKg when the right-handed 
squarks have a mixing. 
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Figure 3: The 

4 = -m,/(47r2)(ll/21)Im[(<5S)23Ci 



correlation 



between 
Here, {5 



ds and S^Ks 
= -0.04 and aTg[C^ 

conies from the matrix element of chromomagnetic moment in B ^ (pKg. The dashed 
(dotted) line is the upperbound on dg from the EDM of ^^^Hg atom (neutron). 
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